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Abstract 

We propose an improved model explaining the occurrence of high stresses due to the 
difference in specific volumes during phase transitions between water and ice. The un- 
knowns of the resulting evolution problem are the absolute temperature, the volume incre- 
ment, and the liquid fraction. The main novelty here consists in including the dependence 
of the specific heat and of the speed of sound upon the phase. These additional nonlinear- 
ities bring new mathematical difficulties which require new estimation techniques based 
on Moser iteration. We establish the existence of a global solution to the corresponding 
initial-boundary value problem, as well as lower and upper bounds for the absolute tem- 
perature. Assuming constant heat conductivity, we also prove uniqueness and continuous 
data dependence of the solution. 

Introduction 

In the present contribution we prove the well-posedness of an initial-boundary value prob- 
lem associated with the following system coupling a quasi-linear parabolic internal energy 
balance (for the absolute temperature 9) with an integro-differential equation for the rel- 
ative volume increment U , and a differential inclusion ruling the evolution of the phase 
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variable x &s follows: 

c(x)ei(0)i-divKx)V0) = d{x)xtUi{0)-eim 

+ uU^-peut + i{e)xl-L^xt, (i.i) 

vUt + X{x){U -a{l-x))- P{0-e,) = gog{xs-Cr)-Kr{Po{t) + Un{t)), (1.2) 
-liO)Xt -^iU-a{l- x)? - aXix){U - a(l - x)) 

G c'ix) Um - Wc)) dl{x) (1.3) 

with Un{t) = U{x, t) dx . In the previous paper |12j we have already given a motivation 
and a complete study of equilibria for this system, which models the water freezing in an 
elastic container, taking into account differences in the specific volume, specific heat and 
speed of sound in the solid and liquid phases. The derivation of the system from physical 
principles and the meaning of the symbols will be explained below in the next Section [2l 
Here, we describe the mathematical difficulties and comment on previous results related 
to this type of systems. 

There is an abundant classical literature on phase transition processes, see e.g. the 
monographs [2], [5]) |14j and the references therein. It seems, however, that only few 
publications take into account different mass densities/specific volumes of the phases. In 
[B], the authors proposed to interpret a phase transition process in terms of a balance 
equation for macroscopic motions, and to include the possibility of voids. Well-posedness 
of an initial-boundary value problem associated with the resulting PDE system is proved 
there and the case of two different densities gi and Q2 for the two substances undergoing 
phase transitions has been pursued in [?]■ 

Here, we deal exclusively with physically measurable quantities. All parameters have 
a clear physical meaning and the derivation is carried out under the assumption that the 
displacements are small. This enables us to state the system in Lagrangian coordinates 
(cf. [7] for a different approach to the subject). 

The present model has been previously studied in [TUj and |llj under the assumption 
that the speed of sound and the specific heat are the same in solid and in liquid. In 
terms of the system (|l.ip - (|1.3p . this corresponds to choose constant functions A(x) = A 
and c(x) = c. For this particular case, we have proved in [10] and [llj the existence and 
uniqueness of global solutions, as well as the convergence of the solutions to equilibria. In 
reality, the specific heat in water is about the double, while the speed of sound in water 
is less than one half of the one in ice. The main goal of this contribution is to give a 
well-posedness result for a boundary value problem associated with ()l.ip - ()1.3p including 
these dependences into the model. The main result is stated in Section [H The dependence 
of speed of sound and of the specific heat on the phase is expressed in terms of additional 
nonlinearities in the equations which have to be suitably handled. Moreover, here we 
also generalize the results of |1U] and [llj allowing for non constant external pressure 
and temperature. Finally, we proceed here with a different technique for the proof of 
existence of solutions with respect to [TO] and [TT] . Since the contraction argument does 
not work in our situation, we discretize in time our problem (cf. Subsection [5Tj, preparing 
thus necessary tools for future numerical investigations on this model, and prove the 
convergence of the scheme. The uniqueness and continuous dependence of solution on the 
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data is proved in Section [6] following the idea already exploited in where we deal with a 
quasi-linear internal energy balance equation coupled with a vectorial and nonlocal phase 
dynamic. The main estimates are obtained here by means of the energy inequality which 
still holds true at the discrete level (cf. Subsection 15. 3|) . Finally, it is worth noting that a 
time dependent positive lower bound for the -component of the solution independent of 
the time step is established on the time discrete approximation in Subsection 15.21 while 
we obtain a uniform in time upper bound on the solution 6 by means of a proper Moser 
estimate (cf. Subsection 15. 4p . 

2 Balance equations 

Referring to [12] for the complete deduction of the model, we consider a liquid substance 
contained in a bounded connected container $7 C M'^ with boundary of class C^'^ . The 
state variables are the absolute temperature > 0, the displacement u G M'^ , and the 
phase variable x £ [0) 1] • The value x = means solid, x = 1 means liquid, x £ (0) 1) is 
a mixture of the two. 

We make the following modeling hypotheses. 

(Al) The displacements are small. Therefore, we state the problem in Lagrangian coordi- 
nates, in which mass conservation is equivalent to the condition of a constant mass 
density go > 0. 

(A2) The substance is isotropic and compressible; the speed of sound and the specific heat 
may depend on the phase x- 

(A3) The evolution is slow, and we neglect shear viscosity and inertia effects. 

(A4) We neglect shear stresses. 

(A5) The liquid phase is the reference state, and the specific volume Vi of the solid phase 
is larger than the specific volume of the liquid phase. 

We thus consider the evolution system 

-diver = fyoi , (2.1) 
Qoct + divq = criSt, (2.2) 

-Jo{0)Xt G dj, (2.3) 

consisting of a the mechanical equilibrium equation (|2.ip . energy conservation law (|2.2p . 
and a phase dynamic equation ()2.3p . where the coefficient 70 determines the speed of the 
phase transition. By (A4), the stress has the form cr = —pd and the scalar quantity 

p := -vet : 6 - \{x){e :6-a{l- x)) + P{0 - Oc) (2.4) 

is the pressure. Here > is a volume viscosity coefficient, A(x) is the Lame constant, 
which may depend on x by virtue of (A2), a = {yi — Vw)/Vw is a positive phase expansion 
coefficient by (A5), while (3 is the thermal expansion coefficient, which is assumed to be 
constant, and iyoi is a given volume force density (the gravity force) 

^vol = -Qog^3 , (2.5) 
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with standard gravity g and vector ^3 = (0, 0, 1) . 

We denote by e the specific internal energy, s is the specific entropy, and q is the 
heat flux vector that we assume for simplicity in the form 

ci=-K{x)Ve (2.6) 

with heat conductivity k{x) > depending possibly on x- 
We assume the specific heat cy(x,6') in the form 

cv{x,9) = co{x)ci{e). (2.7) 

This is still a rough simplification, and further generalizations are desirable. According 
to [3 Chapter VI] or [T31 Section 5] , the purely caloric parts Ccai and Scai of the specific 
internal energy and specific entropy are given by the formulas ecaKX;^) = co{x)ei{9) , 
Scai{x,0) = co{x)si{9), with 

ei{e)= [\iir)dr, = f dr . (2.8) 



Then, the specific free energy / = e — 6s satisfies the conditions = Qodsf , s = 
—dgf. With a prescribed constant latent heat Lq and freezing point > at standard 
atmospheric pressure Pstand , the specific free energy / necessarily has the form 

/ = coix)fiie) + ^{e:S-a{l-x)f (2.9) 
-^{9- e,)e : «5 + Lox f 1 - ^ V Kx) , 

go V c'c/ 

where ^ 

f^{e) = ei{9)-9si{e) = I ci{r)(l--\ dr, (2.10) 



and / is a arbitrary function of x (integration "constant" with respect to 9 and e ) . We 
choose / so as to ensure that the values of x remain in the interval [0, 1] , and that the 
phase transition under standard pressure takes place at temperature 9c ■ More specifically, 
we set 

fix) = LoI{x)-co{x)fi{Oc). 

where / is the indicator function of the interval [0, 1] . 

For specific entropy s and specific internal energy e we obtain 

s = -dgf = co{x)si{9) + ^e:6 + ^x, (2-11) 

£'0 9c 

e = co{x){ei{9)-h{9c)) + ^{e:S-a{l-x)f + -9ce:d + Lo{x + I{x)).{'^-^'i) 

^QO QO 

The equation for the phase x is obtained by assuming that —xt is proportional to 
d^f with proportionality coefficient (relaxation time) 70 (^) > 0, where dy. is the partial 
Clarke subdifferential with respect to x ■ 

Then, the equilibrium equation (|2.ip can be rewritten in the form Vp = f^ol 1 hence, 
as is connected, 

p{x,t) = P{t)-Qogx3, (2.13) 
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where P is a function of time only, which is to be determined. Recall that in the reference 
state £ : S = Et : S = , x = and at standard pressure Pstand i the freezing temperature 
is 9c ■ We thus see from ()2.4p that P{t) is in fact the deviation from the standard 
pressure. We assume also the external pressure in the form P^xt = Pstand + Po with a 
given deviation po{x, t) . The normal force acting on the boundary is {P{t) — QoQ X3—pQ)n, 
where n denotes the unit outward normal vector. We assume an elastic response of the 
boundary, and a heat transfer proportional to the inner and outer temperature difference. 
On dil. , we thus prescribe boundary conditions for u and 6 in the form 

iP{t) - Qogxs-po{x,t))n{x) = k{x)u{x,t) , (2.14) 
(l{x,t)-n{x) = h{x){9 -erix,t)) (2.15) 

with a given symmetric positive definite matrix k(x) (elasticity of the boundary), positive 
functions h{x) (heat transfer coefficient), and 9r{x,t) > (external temperature). This 
enables us to find an explicit relation between divu and P. Indeed, on dfl we have by 
(|2.14|) that u • n = {P(t) — gog x^ — pq{x , t))\i.~^ (x)n(x) • n(x) . Assuming that k~^n • n 
belongs to L^{dQ,), we set 



— = f k ^{x)n{x) ■ n{x) da{x) , = Kr j k ^{x)iy{x) ■ ^.{x) X'ida{x) , (2. 
Jdo. Jan 



Kr 

and obtain by Gauss' Theorem that 
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Unit) := J^diyu{x,t)dx = ^iP{t) - gogCr) - Poit) , (2.17) 

where Po{t) = jQ^pQ{x,t)]<.~^{x)n{x) • n{x) da{x) . Under the small strain hypothesis, 
the function divu describes the local relative volume increment. Hence, Eq. (|2.17|) estab- 
lishes a linear relation between the total relative volume increment Unit) and the relative 
pressure P{t) —po{x,t). We have e : S = divu, and thus the mechanical equilibrium 
equation (f2l^ . due to ([23]) and (f2l71) . reads 

i/div ut + A(x)(div u - a(l - x)) -^{9-9^) + QogiCr - X3) = -Kr{Po{t) + Un{t)) . (2.18) 

As a consequence of ()2.6p . ()2.9p . and ()2.12p . the energy balance and the phase relax- 
ation equation in (|2.2p - (|2.3p have the form 

QoCo{x)ei{9)t - div(K(x)V0) + 0oc'o{x)Xt{ei{9) - h{9)) 

Q 

= uidivutf - P9dwut + QoJo{9)xt - QoLo^Xt, (2.19) 



-Qo7o{0)xt - ^^^^^(divu - a(l - x)f - aA(x)(div u - a(l - x)) 

G Qoc'Ax) ifiiO) - h{9c)) + 0OLO (1 - ^) + 91 ix) ■ (2.20) 
For simplicity, we now set 

C/:=divu, c{x) := Qocoix) , l{0) := Qo7o{9) , L := qqLq . (2.21) 
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Note that mathematically, the sub differential dl{x) is the same as g()LQdI{x) ■ The 
system thus reduces to the system (|l.ip ~ (|1.3p of three scalar equations - one PDE and two 
"ODEs" for three unknown functions 9,Xi ^^(^ U , with boundary condition ()2.15p . ()2.6p . 
Assuming that a solution to (|l.ip - (|1.3p is known with U G L'^{fl x (0,r)), we find the 
vector function u by defining first <I> to be the solution to the Poisson equation A<1> = U 
with the Neumann boundary condition V<l>-n = {KrUn{t)+ gog{Cr—X3))'k.~^ {x)n{x)-n{x) . 
With this <I> , we find u as a solution to the problem 

divu = in nx (0, T) , (2.22) 

u • n = 



(u + V^> - (KrUn + Qog{Cr - X3))k-^n) x n = 



on dn X (0,T) , (2.23) 



and set u = u + V<I>. Then u satisfies a.e. in Q the equation divu = U , together with 
the boundary condition ()2.14p . that is, u = {KrUn + Qog{Cr — xs))'k^^n on dO,. 

For the solution to (|2.22P " (j2.23p . we refer to [H Lemma 2.2] which states that for each 
g G H^^'^{dil.)^ satisfying Jq^ g • nda{x) = there exists a function u G H^{Q)^ , unique 
up to an additive function v from the set V of divergence-free H^{Q) functions vanishing 
on dil., such that divu = in ft, u = g on dil.. In terms of the system (|2.22p - (|2.23p . 
it suffices to set g = ((V<I> — (KrUfi + Qog{Cr — a^3))k~"'^n) x n) x n and use the identity 
(b X n) X n = (b • n) n — b for every vector b. Moreover, the estimate 

mf ||u + v||^i(f^) < C ||g||Hi/2(9n) < C\\<^\\H2^n) (2.24) 

holds with some constants C, C . The required regularity is available here by virtue of 
the assumption that Q is of class C^'^, provided belongs to H^/'^{dVl) . Note that a 
weaker formulation of problem ()2.22p - ()2.23p can be found in [H Section 4]. 

Due to our hypotheses (A3), (A4), we thus lose any control on possible volume 
preserving turbulences v G F . This, however, has no influence on the system ()l.ip - ()1.3p . 
which is the subject of our interest here. Inequality (|2.24p shows that if U is small in 
agreement with hypothesis (Al), then also v can be chosen in such a way that hypothesis 
(Al), interpreted in terms of , is not violated. 

3 Energy and entropy 

In terms of the new variables 6,U,x, the densities qqb, qqs of energy and entropy can be 
written as 

Qoe = c{x)iei{e)-hm + ^{U-a{l-x)? + /30cU + L{x + Iix)), (3.1) 
Qos = c{x)si{e) + ^x + PU. (3.2) 
The energy functional has to be supplemented with the boundary energy term 

m) = ^ (unit) + Po{t) + , (3.3) 
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as well as with the gravity potential —Qogx^U . The energy and entropy balance equations 
now read 



Qo{e{x,t) - gx3U)dx + Erit) 



h{x){er{x,t) -e)da{x) 



+ Kr{Po)t{t) (Un{t) + Po{t) + , (3.4) 

Qost + div- = ^2 + —0-Xt+-gUt > 0, (3.5) 

^ [ g,s{x,t)dx = [ ^{er{x,t)-9)da{x) (3.6) 
di Jn Jdn 

The entropy balance ()3.5p says that the entropy production on the right hand side is 
nonnegative in agreement with the second principle of thermodynamics. The system is 
not closed, and the energy supply or the energy loss through the boundary is given by the 
right hand side of ()3.4p . 

We prescribe the initial conditions 

e{x,0) = 9^{x) (3.7) 
U{x,0) = C/°(x) (3.8) 
X(x,0) = x°(x) (3.9) 

for x E and compute from (|3.ip - (|3.2p the corresponding initial values , E'p , and 
for specific energy, boundary energy, and entropy, respectively. Let and 5'^ denote, 
respectively, E^ = J^goe^{x)dx, = qos^{x) dx . From the energy end entropy 
balance equations (j3.4p . ()3.6p . we derive the following crucial (formal for the moment) 
balance equation for the "extended" energy QQ{e — Oys), 9r being a suitable positive 
constant: 

c(x)(ei(0) - fiiOc)) + ^{U - ail - x))') ix,t)dx 

+ / i/SOcU + Lx- gogx3U)ix,t)dx 
Jn 

+ ^(unit) + Poit)+''''^'^' 



+ /* / ^iO-0rix,m0-Or)daix)d^ 
Jo Jan 

E^ + E^- BrS^ + Or [cix)si (0) + jX + (iU^ ix. t) dx 

+ j'^KriPo)tii){unii) + Poii) + ^^ d^ (3.10) 



We assume that both c{x) and A(x) are bounded from above and from below by positive 
constants. The growth of 51(6*) is dominated by 61(6*) as a consequence of the inequahty 

Sl{e)-Siie*) 1 

77T^ TTTT ^ TT > & > . 

ei{e) - ei{e*) - e* 

We win use the relation ()3.10p to get an upper bound for the solution on the whole time 
interval (0, oo) . From the identity 

he -a){e-h) = \{9- V^f -{Vb- (3.11) 

U a 

for all ^, a, 6 > 0, it follows that we find a constant C > independent of t such that for 
all t > we have 

I (em + u^) *) d. + /' I + ^ + ^) (-« ^^df 

+ I' L^i"- V^tSm) d<r(l-) d{ < C , (3.12) 

provided we assume that 

/ / h{x) (y^9r{x,t) - y%) da{x)dt<oo, / \{Po)t{t)\dt < oo . (3.13) 
Jo JdQ ^ ^ Jo 

4 Main results 

We construct the solution of ()l.ip - ()1.3p by a combined truncation and time discretization 
scheme. The method of proof is independent of the actual values of the material constants, 
hence we choose for simplicity 

L = 2, ec = a = /3 = i^ = Qo = l. (4.1) 

We consider the following assumptions on the data. 

Hypothesis 4.1. Assume that there exist positive constants , c* , c, c, A, A, k*, A*, 
7* such that 

(i) c convex, c G C^'H[0, 1]) , < < c{z) , < c < c'{z) < c, for aU z G [0, 1] ; 

(ii) ci G C70(R+), ci{e) > c* for 9 > I, limg^^ci{e)/9 = oo, ei(^) := f^ci{r)dr, 
Jo ci{r)/rdr < oo , ci (r) /r^ dr = oo ; 

(iii) A convex, A G C^'^i[0, 1]) , < A < X{z) < A, > A'(z) > -A* for aU z G [0, 1] ; 

(iv) K, G C^'HIO, 1]) , < < k(z) for aJi z G [0, 1] ; 

(v) /i G L°°(9r2) is a non-negative function; 

(vi) 7 G C°'HlR"^) , < 7* < 7(r) for aii r G IR+ . 

The liquid phase does not persist for very large temperatures and the behavior of 
ci{9) as — 7- oo thus cannot be experimentally verified. We nevertheless believe that the 
growth condition (ii) in Hvpothesis 14. 1 1 is not completely meaningless taking into account 
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the fact that in the interval between 273 and 373 K (0-100° C), the function ci{9) is 
convex with a minimum at 35 ° C 0- 
We introduce the following notation: 

A{U, X, X, t) := X{x){U - 1 + x) + Kr{Un{t) + Po{t)) + ^(Cr - ^3) + 1 , (4.2) 
B{X,0) := c{x){fi{0) - MOc)) - 29, (4.3) 

m x):=^iU-l + xf + KX){U - 1 + x) + 2 . (4.4) 
System (|l.ip - (|1.3|) with boundary condition (|2.15|) then can be written in the form 

c{x)ei{6)tw{x) dx + / k(x)V0 • V7i;(x) dx = / h{x){6r{x,t) — 9)w{x) da{x) 
Jn Jan 

UtA{U, X, X, t) + xt{C{U, x) + c'{x){ei{9) - h{e,))))w{x) dx , (4.5) 

Ut-e = -A{U,x,x,t), (4.6) 
l{0)xt + B{x, 9) + dl{x) 3 -C{U, x) , (4.7) 

where (|4.5p is to be satisfied for all test functions w € VF"'^'^(0) and a.e. t > 0, while (|4.6|) - 
()4.7p are supposed to hold a.e. in the space-time cylinder that we denote VLt := O x (0, T) 
for r > 0, Ooo := X (0,00). 

In this section we prove the following existence and uniqueness result. 

Theorem 4.2. Let Hypothesis be satisfied, and let 6r G H^{0,T] L^{dQ)) such that 
< 6^ < 9r < 9* , and Pq G 1^-'^'^(0,T) be given functions. Let the initial conditions in 
(ISZD-dSI) be such that 

9° eW^'^{n)nL^{n), 0<9^<9^{x) <9* a.e., 
C/°,X° e T^^'2(f])nL°°(0), 0<x°(s)<l a.e. 

Then there exists at least a solution [9,U,x) to ^^-^T^, 71 )-( fXP)) . and constants 
9\T) > 9\T) > O such that 

9\T) < 9{x, t) < 9\T) for a.e. {x, t) £ Qt , (4.8) 

X G [0,1] a.e., U,Ut,Xt G L^i^r), 9t G L^i^r), VU,VxyO G L-(0, T; L2(f])) . // 
moreover condition ()3.13p is satisfied, then the solution exists globally, and 9'^{T) can he 
chosen independently of T . Finally, if k{9) = R £ is constant, then the solution is 
unique, and its L^ -norm depends continuously on the data. 

Remark 4.3. Let us note that we could prove our existence result assuming that k = 
k{9,x) = ki{9)k2{x) with the same techniques. Moreover, also uniqueness would hold 
true in case k = ki{9) with an appropriate modification of the boundary condition by 
means of the standard Kirchhoff transformation technique. 



see 



http : //www . engineeringtoolbox . coin/water-thermal-properties-d_162 . html 
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5 Existence proof 



We proceed as follows: We truncate from above the functions depending on 6 in ()4.5p - 
()4.7p . and discretize the system in time. For the discrete system, we derive upper and 
lower bounds that enable us to let the time step tend to and prove the existence of a 
solution to the truncated problem. Finally, we prove a time dependent lower bound and 
a uniform (in time and w.r.t. the truncation parameters) upper bound on 6 , so that the 
truncation can be removed, and this will conclude the proof of existence of solutions. 

5.1 Approximation and discrete energy estimate 

We introduce, for ^ G M, i? > 0, the functions 

QR{e)=mln{e^,B{R)], B{R) = R^I\nAn{ei{R),\h{R)\]f'\ (5-1) 
cf{e) = ci{QR{9)), (5.2) 

ef (^) = / cf (r) dr, (5.3) 







.fW^/^d. (5.4) 







fhO) = ef (^) - QR{e)sm = J' cf (r) (l - ^) dr. (5.5) 

In the rest of the proof the following relations, which directly follow from the above 
definitions and from Hvpo. Hm (ii). play an important role: 

If 9 < B{R) then ef{9) = ei(^), sf{9) = si{e), f^{9) = MB), (5.6) 
U9> B{R) then ef (0) = ei{B{R)) + ci{B{R)){9 - B{R)), 

sf{e) = s^{B{R)) + ^^c,{B{R)){9 - B{R)), ff{9) = h{B{R)), (5.7) 
If ^ > i? then ef{9) > ei{R) > 0, /f (0) < fi{R) < 0, (5.8) 

Hm ^ = - hm ^ = - hm ^ = -co, (5.10) 

\h{R)\ ei{R) B{R) 

1™ no/m = lim 75977^=00) „ = oo. (5.11) 

We now introduce the time-discrete version of ()4.5p - ()4.7p . For an arbitrary n G N, 
we define the time step r = T/n. Choosing a constant cr G M.^ depending on R, which 



10 



we specify below, we look for a solution {{Of^, Uk,Xk)}k=i scheme 

- / c{xk){ef {ek)-ef{ek-i))wix)dx+ [ K(xfc-i)V0fc • Vu;(x) dx 



r Jq 



CRiOkOt - Ok-iOt^i)w{x)dx + / h{x){ek-9kr)w{x)da{x) 

an 



c'iXk) Xfc(ef - fi{6c))w{x) dx 

(UkAk{Uk, Xk, Xfc-i, x)+ XkCk{Uk, Xk, Xk-ifj w{x) dx , (5.12) 

Uk - QR{Ok~i) = -■Ak{Uk,Xk,Xk-i,x) , (5.13) 
l{Ok~i)Xk + Bk{Xk,Ok^i) + dI{xk) 3 -Cfc(C/fc,Xfc,Xfc-i), (5.14) 

where 



Ak{Uk,Xk,Xk~i,x) := X{xk-i)iUk - 1 + Xfe) + Kr{Ukn + Pk) + giCr - X3) + 1 , (5.15) 
BkiXk, 0k-i) ■■= c'(xfc)(/f (^fe-i) - fi{Oc)) - 2QR{6k-i) , (5.16) 

Ck(.Uk,Xk,Xk-i) ■■= ^^^{Uk-l + Xkf + X(.Xk-i)(.Uk-l + Xk) + 2, (5.17) 



with 9kr = Ori-, kr), pk = Poikr) , Ukn = Un{kT), and with initial condition {9o,Uo,xo) 

= {6^,U^,x^)- The symbol 'f^ denotes the discrete time derivative {fk — ^k-i)/T for a 
generic sequence {^Pk} ■ Eq. (|5.12p has to be satisfied for all test functions w G W^''^{Q,), 
while (|5.13|) - (j5.14p are supposed to hold a.e. in . 

It is easy to see then that the latter semi-implicit scheme has a unique solution. Indeed, 
at each time step, we assume that 9^-1, Uk-i,Xk-i are known, and find Uk,Xk satisfying 
(|5.13p - (|5.14p . For r sufficiently small, (|5.13p - (|5.14p is an algebraic system for {Uk,Xk) of 
the form ^{Uk,Xk) = ^fc with $ : — >• strictly maximal monotone, hence it admits 
a unique solution. Finally, we insert Uk and Xk in ()5.12p and solve the resulting coercive 
elliptic equation, obtaining in that way the desired solution {9k,Uk,Xk) ■ 

Then, we note that the total energy balance still holds true for the discrete system. 
Indeed, we take (|5.12p with w = 1, and denote = ef (0^) — fi{9c) ■ We have 



- / cixk){e?{9k) - efi9k-i)) + ixk - Xk-i)cixk)ie^i9k-i) - fii9c))dx 
T Jq 

= - c{xk){Ek - Ek-i) + ixk - Xk-i)c{xk)Ek-i dx 
T Jn 

= - c{xk)Ek - (c(xfc) - iXk - Xk-i)c' {xk))Ek-i dx , 



and using the fact that E^ > and that c is convex (cf. Hvpo. HTT] (i)). we get 

Z / <Xk){e^{9k)-e^{9k.^)) + {xk-Xk-i)c'{xk){ef{9k-i)-h{9c))dx 

> - [ {c{xk)Ek - c{xk-i)Ek-i) dx. (5.18) 
T' Jn 



r Jq 
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Set now Sk = Uk + Xk ~ ^ ■ Then, we obtain 

iXk+ Uk)KXk-i){Uk + Xk- 1)+ Xk^^Y^iUk + Xk- I? dx 

If 1 

= - / ^^(Sfe - 5fc_i)A(xfe-i) + 7T'S'|(Xfc - Xfc-i)A'(xfc)dx 

- h X^^' " SLi)^{Xk^i) + SUxk - Xk~i)X'{Xk) dx 

If 1 
= ^ J^Sl (A(xfc-i) + A'(xfc)(xfc - Xk~i)) - ^5'Li-^(Xfe-i) dx . 

Using now the convexity of A (cf. Hypo. UT] (hi)), we get 

/ iXk+ lJk)KXk^i){Uk + Xk- 1)+ Xk^^^^iUk + Xk- I? dx 
Ju ^ 

>^ f SlXixk) - ^LiA(xfc-i) dx . (5.19) 



2r „ 

Hence, from (|5.12p . using (jS.lSp and ()5.19p . we obtain 

i ^ (^c{xk)Ek - c{xk-i)Ek^i + \ {SlXixk) - Sl^.Xixk^i))^ dx 

+ Ukn{Kr{Ukn+Pk) + gCr) - / gx^Ukdx 

Jq 

+ I (Uk + 2Xk + CRi9ke+ - ek-i9+_^)) dx+ [ h{x) {Ok - Okv) da{x) 
Jn ^ ' Jan 

<0. (5.20) 
Summing now (|5.20p over k = 1, . . . ,m , 1 < m < n, we get 

j (^c{Xm)Em + ^5'^A(Xm) " gXsUm + Um + Xm + CRTdmO£^ dx 

+ ^ + +Tf2j^^h{x){ek-ekr)da{x) 



< (c(xo)^o + \sIX{xq) - gxsUo + C/o + Xo + cr^o^^) dx 
+ ^(uon+Po + ^] + Kr^^lpk - Pk-i\ max 

Z \ Ap / ^ — ^ 0<fc<»n 



Ukn +Pk + ^ 
Ar 



(5.21) 



Using the fact that r /an ^(a^)^fcr(x) da(x) < C{T), El=i\Pk-Pk-i\ < C(r), with 
a constant C(T) independent of r and R, we check that the left hand side of (j5.2ip is 
bounded independently of r and R. Consequently, all terms in Eq. (|5.13p are bounded 

by a multiple of (1 + B{R)) . Similarly, multiplying (|5.14p by Xk and using the fact that 

Xk^k > for all ^k G dl{xk) , we obtain the estimates 



\Uk\ + \Uk\ < CiT){l + B{R)) 

\Xk\ < C{T){l + B{R) + B\R) + \h{B{Rm 



a.e. (5.22) 
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5.2 Lower bound for 6i 



Here we derive a lower bound for the approximated absolute temperature 9k ■ We first 
rewrite SUTHh for w £ W^'^{n) , w>0 a.e., using (j5.13p - (|5.14p . in the form 

- / c{xk){e^i0k)-e^{ek-i))w{x)dx+ [ K(xfc_i)V4 • Vw;(x) dx 
T Jn Jn 

+ / {cR{9k9+ - 9k-i9+^^)) w{x) dx+ [ h{x){9k - 9kT)w{x) da{x) 
Jn JdQ 



> / Uk{Uk-QRi9k-i))wix)dx 

+ [ Xkil{0k^i)Xk + c'{xk){ff{0k~i)-ef{9k-.i))-2QR{9k^^))w{x)dx, (5.23) 
Jn 

where we have used again the fact that 

- XkC{Uk,Xk,Xk-i) > Xk{-f{9k-i) Xk + 13{xk,0k-i)) 

by definition of the subdifferential. The right hand side of (j5.23p is bounded from below 
by a negative multiple (depending on R) of ^^-i^fcLi- We can now choose cr in ()5.12p 
sufficiently large in order to get the following inequality for all w G W^''^{Q,) , w > a.e.: 



1 

r Jn 



c{Xk) {enOk) - ef{9k-i)) w{x) dx + / K{xk~i)V9k ■ Vw{x) dx 



+ / h{x){9k - 9kv)w{x)da{x) > -cr I 9k9'^w{x)dx. 
Jan 



(5.24) 



We now compare this inequality with the constant decreasing sequence {vk} defined 
recurrently as 

-c* (ef (ufc) - ef{vk-i)) = -CRvl, vq := 9^ . (5.25) 
We write ()5.25p . adding the zero term — div (A;(xfc-i)Vi;fc) , in the form 



(5.26) 



Subtracting ()5.24p from ()5.26p and testing the difference hy w = Hf,{vk — 9k), where Hi; 
is the regularization of the Heaviside function H , 



-J {ei (vk) - ei {vk-i)) w{x) dx + j Hi{xk-i)^Vk ■Vw{x) dx 
= —cr I v\w{x)dx. 



He{v) 



if t; < 
v/e if u G (0, e) 

1 if t; > e 



(5.27) 



we obtain, since Vk < Vk-i , that 

/ c{xk) {{efivk) - ef{vk-i)) - {ef{9k)-ef{9k^i)))H,{vk-9k)dx<0. (5.28) 
Jn 
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Assume that O^-i > v^-i (this is true for k = 1). For e\0, (|5.28|) yields 6k>Vk, and 
by induction we get 9k > > Vn for all k = 1, . . . ,n. By (|5.25|) . we have 

ei(wfc) - ei{vk-i) = -CtvI 

with C = cpt/c^: . Under Hypo. 14.11 (ii), we have that the function G{z) = — JJ" ^^^^ ds = 
+00. Then, G is increasing in (0, f o] , G(0+) = — oo, G{vo) = 0. Moreover, by the Mean 
Value Theorem, there exists Sk G [vk,Vk-i] such that 

Gjvk) - Gjvk.,) _ G'jsk) _ 1 ^ 1 
ei{vk) - ei{vk-i) ci(sfc) 4 ~ ^fc ' 

hence G{vk-i) — G{vk) < Gt, that is, G{vn) > —Gnr and so 9k > Vn > G~^ (—Gtit) = 
G~^{—GT) =: 9^{T) . This concludes the proof of the lower bound for 9k ■ 



5.3 Estimates 

Now, we perform the estimates we need in order to pass to the limit as r \j in (j5.12p - 
()5.14p . The right hand side of ()5.12p is bounded from above, by virtue of (|5.22p . by 
C(T, i?)(0fc_i + l) , where C(T, R) is, here and in the sequel, any sufficiently large constant 
depending only on T and i?, and independent of k and r. Testing (|5.12p hy w = 
9k — 9k-i , we obtain 

^ (^clxfc) (ef - ef (^fc-i)) i9k - ^fc-i) + f^ixk^i)V9kV{9k - dx 

+ I CR{9k - 9k^i)\9k + 9k^i) dx+ [ h{x){9k - 9kv){9k - 9k-i) da{x) 
Jn Jan 

<G{T,R) [ \9k-9k-i\{9k-i + l)dx. 
Jn 

Using the lower bound for 9k , and choosing x-i = Xo > we get 

-p^, [ \9k - 9k^i\^ dx + [ (K(xfc-i)|V0,.|2-K(xfc_2)|V^fc_i|2) dx 
T^iK-L ) Jn Jn 

+ [ h{x) {{9k - 9kv? - {9k-i - 9i^k-i)T?) da{x) 
Jan 

< (K(xfc-i) - K(xfc-2)) |V6ifc_i|^ dx+ / h{x)\9kr-9(^k-i)r\\9k-9kr\da{x) 
Jn Jan 

+ tG{T, R) [ {9k^i + if dx 
Jn 

< tG{T, R) (1 + + |V0fc_i|2) dx + h{x) (l ferl' + {9k - 9kr?) da{x)^ , 



where Gi is a positive constant depending on T but not on r . The elementary inequality 
9k — 9l_i < ^{9k — 9k-if + ■^{9k + 9k-if enables us to rewrite the above inequality in 
the form 

Qk - Qk-i < TG{qk + Qk-i + bk-i) , (5.29) 



14 



with 



hk = l+ / h{x)\dkT? da{x) , 
J an 

and C := C{T,R). Inequality (|5.'29p is equivalent to 

1 + tC tC 

qk < ^ _ + _ ^(j ^k-i , (5.30) 

which yields 

k ^ k-l , .k-l-j ( fc-i 



Qk 



< 



T^TcJ '^o + r J:^,^ j < e--Uo + r5:6, (5.31) 

^ j=0 ^ \ j=0 j 

holding true for r < 1/(3C) . We conclude for all m = 1, . . . , n that 

/ |0fe-0fc_ipdx+ / |V0„,pdx+ / /i(x)(0™-0^r)' da(x) <C(r,i?). (5.32) 

Then, we introduce the piecewise constant and piecewise linear interpolants, for t S 
[{k — 1)t, kr) , k = 1, . . . ,n, hy the formula 

^M(x,t) = ^fe_i(x), 9^^\x,t) = ek{x), e^^\x,t) = 0k-i{x) + {t-{k-l)T) hix), (5.33) 
with a similar notation for U , Xi , and Pq . In particular, we set 

eW(x, t) = ef (0fc_i(x)) + i(t - (A; - l)T)(ef (^^(x)) - ef (0fe_i(x))) . 

r 

The previous estimates give immediately that 

§1''^ bounded in L'^{0,T; L^{n)) , 
V^W bounded in L°°{0,T; L^{n)) , 

£ j {\&'^ - ^^"^ + 1^^"^ - ^^"^1^) {x,t) dxdt < C{T,Ry . 
Letting r tend to and passing to subsequences if necessary, we get the convergences 

(5.34) 



^ e strongly in C°([0,r];L2(17)), 

e^''^ 0,6'^^') 9 strongly in L2(0,r;L2(O)), 



e'^^'' Ot weakly in L2(0,T;L2(J7)), 

V^(^) ^ weakly* in L°°(0,r;L2(J])) 



Now we estimate Vxfc and Vf/^ as follows. Prom Eq. ()5.13p it follows that 

{Uk{x)- Uk{y)){Uk{x) - Uk{y)) < C{T,R){\Uk{x) - Uk{y)\' + \xk{x) - Xk{y)\^ 
+ \xk-i{x) - Xfc-i(y)P + l^fc-i(^) - 0k-i{y)\'^ + \x- yp) 
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and, analogously, from (|5.14|) . we obtain 

Xk{y)){xk{x) - xk{y)) < C{T, R)[\Uk{x) - Uk{y)\^ + \xk{x) - Xk{y)\^ 
+ \xk-i{x) - Xfc-i(y)P + \Ok-iix) - ^fc-i(y)p) . 

Summing up the two previous inequalities, we get 

{Uk{x) - Uk{y)? + {Xk{x) - Xk{y)? < (Uk^iix) - Uk^Mf + {Xk~i{x) - Xk^M? 
+ TC{T,R)[{Uk{x) - Uk{y)f + {xk{x) - xk{y)? + {Uk~i{x) - Uk-i{y)f 
+ {Xk-i{x) - Xk-i{y)? + {ek-i{x) - 9k-i{y)f + {x- yf) . 

We are again in the situation of Eq. ()5.29p . with = {Uk{x) — Uk{y))'^ + {Xk{x) — Xk{y))'^ , 
hk = {9k~i{x) - 9k^i{y)f + {x- yf. Hence, by (I5.3ip. we obtain 

{Uk{x)-Uk{y)? + {Xk{x)-Xk{y)? 

k-l 

< C(r,i?)((C/o(x) - Uo{y)f + ixoix) - xo{y)f + TY,iOj{x) - e,{y)f + {x- yf) . 

j=Q 

Using now the previous estimate on V^"^ , we get 

V[/(^), Vx^^^ bounded in L°°{0,T; L'^{n)) . 
We already know that ul:^\x['^^ are bounded in L°°{Qt)- Furthermore, 

/"^ /" |C7W -^7M|2(x,t) dxdt < C{T,R)t\ 
Jo Jn 

I ^(lx^"^-X^"^l' + lx^"^-x("f)(^,i)dxdt < C{T,Ry, 

so that the convergences ()5.34p take place also for U and x ■ We now rewrite ()5.12p - (j5.14p 

in terms of the functions 6^''^ ,e^^\ e(^) , x^^^ , x^^^ , X^""^ , U^'^'^ , f^^""^ , , H''^ ■ The above 
estimates allow us to pass to the limit as r \ and obtain a solution for the following 
truncated problem 

/ c{x)ef{9)tw{x)dx + / k{x)V9 ■ Vw{x) dx = / h{x){9r{x,t) - e)w{x) da{x) 
Jn Jn Jdn 

- [ {UtA{U,x,x,t)+xt{C{U,x) + c'{x){ef{9)- fi{ec))))w{x)dx, (5.35) 
Jn 

Ut-QR{9) = -A{U,x,x,t), (5.36) 

l{9)xt + Br{x, 0) + dl{x) 3 -C(f/, x) , (5.37) 

where 

Br{x,9) := c'ix){f?{9) - h{9c)) " '^Qr{0) , (5-38) 

A and C are defined in (|4.2p . (j4.4p . and (|5.35p is to be satisfied for all test functions 
w e T^i'2(f]) and a.e. t G (0,r), while (fOBjl - ffOTjl hold a.e. in fir- 
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The next step consists in proving that 9 remains uniformly bounded also from above 
independently of i?, so that the truncation does not become active if R is sufficiently large. 
The argument is based on the following counterpart of the extended energy balance ()3.10p , 



c(x)(ef (^) - /i(0c)) + ^{U - 1 + x)') (x, t) dx 
+ I {U + 2x-gxsU) {x,t)dx + ^(^Un{t)+Po{t) + ^^ 

+ /V ^A0-er{x,mQB{e)-6r)da{x)d^ 
Jo Jan Qr\P) 

= E^ + E^- erS° + er [ (c(x)sf (^) + 2x + U) {x, t) dx 

Jn 

+ 1^ KriPoW) (uniC) + PoiC) + dC , (5.39) 

which holds for every solution to (j5.35|) - (|5.37|) and every t G (0, T) . 

5.4 Uniform upper bound for 9 

We choose R large enough such that B{R) > 6** > . Then {d-Orix, 0){Qr{^) - ^r) > 
{Qr{0) — Oy{x,C)){Qr{(^) ~ 0r) ■ We may therefore argue as at the end of Section [3] and 
obtain from (fO^ for all t G (0, T) that 

< c(i+ / / ;<io(\/«R^7e)-\/ft^)'d<T(i)d{+ / i(F„),({)|d{). (5.40) 

In order to perform the Moser iteration scheme on 9 as in Prop. 3.6], we need first 
to estimate U and Ut in terms of 9 . Rewriting (|5.36p as 

Ut + X{x)U = QR{9) + G{x,t) 

where, by virtue of ()5.40p . G{x, t) is bounded above by a positive constant Go . Denoting 
A(x,t) := JqX(x(x,s)) ds, we obtain the formula 

i7(x,t) = e-^("'*)f/o(x)+ re^("'^)-^("'*)(Qjj(0) + G)(x,e)de. 

JO 

Using Hypo. 14.11 (iii) , we get the estimate 

\U{x,t)\ <\Uo{x)\+ 1^ e^^'-^^Qni9){x,0dC + ^ (5.41) 
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and 

\Ut{x,t)\<\\U{x,t)\+Qnmx,t) + Go. (5.42) 

Now we are ready in order to start the Moser iteration scheme. Choose in ()5.35p 
w^x) = , u = i1^r{9) := {Qr{9) — R)~^, with any p > 1 and with R larger than the 
constants in Hypothesis 14.11 Then 



X){efmtu'' dx + -^^ [ K{x)\Vu'^\^dx+ [ h{x)uP+Us{x) 

<- [ {UtA{U,x,x,t)-xt{C{U,x)-c'{x){ef{e)-Me,))))uPdx. (5.43) 
Jn 

Put Ep{6) = Ci'{r)ilj^{r) dr . Then, we can rewrite (|5.43p as 

(c(x)i?p^(0)), dx + ^^-2 / Kix)\Vu-^\''dx+ [ h{x)uP+Us{x) 
^ * {P+^Y Jn Jan 

< - [ UtA{U, X, X, t)uP dx 



n 

Xt {CiU,x)u^ + c'(x)((ef (0) - Mec))uP - E^{e))) dx . (5.44) 



We now prove that the last integral in ()5.44p is non-positive if R is sufficiently large. First 
of all let us note that, if xt = then it vanishes. Hence, let us consider the case Xi 7^ 0- 
Then, from (fOTD it follows that xt{l{0)xt + Br{x,0) + C{U,x)) = 0, hence 

Xt = -^^{BR{x,e) +C{U,x)) . 

The last integral in (j5.44p is of the form — t^-^i ^ l2dx, where 

h := xt = -c{u,x) + '2QR{9) + d{x){h{ec)- fho)), 

h := C{U,x)u^ + c'{x){{e'^{e)-h{e,))uP-E^{e)) 

= {C{U,x) - c'{x)fi{Oc))nP + c'ix) (ef (0)nf - Ej^{e)) . 

We can now estimate from below the last term as follows 

e^{e)uP - E^{e)= / pef (r)nP-Mr > ei(i?)nP. 
Jo 

We have /2 = if < R, while for 9 > R we have by Hvpothesis I4.1l (i) 

h > -\CiU,x)\+2R + c\ff{R)-M9,)\, 
h > uP{-\C{U,x)\+ch{9,) + cef{R)). 

By virtue of (15:22]) . we have \C{U,x)\ < C{B'^{R) + 1). Referring to (l5TT]) . we conclude 
that there exists Rq> \ larger than all constants in Hvpothesis 14. II such that for R> Rq 
we have in (|5.44p 

Xt {C{U, x)yP + c'(x)((ef (^) - h{9c))uP - E^{9))) dx<Q. 
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Let us fix now R > Rq and continue the Moser estimate, rewriting (j5.44p as follows 
^ (c(x)£^p^(0)), dx + ^ Kix)\W-^ |2 dx + hix)uP+' dsix) 

< - [ UtA{U,x,x,t)uPdx. (5.45) 
Jn 

We have c{x)Ep{u) > ^^u^~^^ , k{x) > . Integrating ()5.45p in time, we obtain, using 
Hypo.[3?T](i),(ii),(iv), as well as the estimates ()5.4ip - (|5.42p and the fact that u{x,0) = 0, 
that 

C*C* f 50+1/ ,N , , 4pK^ 



< 



^f+i(x,i)dx + -^^ / / Vn^ (x,e)dxde 
P + 1 7f7 (P+ 1)^ Jo Jn 

f [ uP{x,o(l + QR{eix,0)+ [\~^^^~'^^QRi9{x,r]))d7]) r{x,Odxd^ . 
Jo Jn \ Jo J 



The function r{x,t) = CA(JJ,x,x,t) , where C is a suitable constant, has norm in 
L°°(0,r;L2(O)) bounded independently of R by virtue of dOOjl . Note that Qnie) < 
u + R. Hence, the function v := u/R satisfies for all p > 1 the inequality 



v' (x, t) dx + 







/o Jn 





2 

(x, ^) dx d^ 



< 



/ / t;*'(x,e) fl + t^(2;,0 + /'^e-^(«-^)7;(x,7?)dr/) r(x,Odxde. 
Jo Jn \ Jo J 



The argument of [TOl Prop. 4.5] yields ||7;||j;^oo(f7y') < C with a constant C independent 
of R and T . Consequently, 

\\QRm\Loo(^nr) < {l + C)R. 

Choosing R sufficiently large such that B(R) > (1 + C)R, we can remove the truncation 
from (j5.35p - (|5.37p . concluding in this way the proof of existence of a bounded solution to 
()4.5p - ()4.7p . If moreover (|8T3]1 holds, then r £ L'^{0,oo; L^{n)) , and the upper bound 
holds globally in iloo • Indeed, the lower bound for 6 in Subsection 15.21 is independent of 
the time step r and is preserved when r \ . 

6 Uniqueness and continuous data dependence 

In this Section, we prove uniqueness and continuous data dependence of solutions under 
the more restrictive assumption that K{r) = S € R"^ for all r € M"*" . 

In what follows, we denote by Rq, Ri, R2, ■ ■ ■ suitable constants that possibly depend 
on T, but not on the solutions. We first rewrite Eq. ()4.5p in the form 



[ {c{x){ei{9)- fi{ec)))tw{x)dx+ [ RV9-Vw{x)dx+ [ h{x){e - er)w{x) da{x) 
Jn Jn Jon 

{UtA{U, X, X, t) + xtC{U, x))w{x) dx , (6.1) 



and denote 6 = 0i- 62, X = Xi- X2, Xo = Xoi - X02 , Oo = doi - do2 , ^^r = On - 0r2 ■ 
U = Ui-U2, Uo = Uoi- U02, e{x,t) = /o*^(x,T)dT, Gr{x,t) = J^^ erix,T) dr , Pq = 
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Poi - Po2- Within the range e\T) < 6* < 6'«(T) and x G [0, 1], \xt\ < C of admissible 
values for the solutions, and, thanks to Hypo.HTTJ all nonlinearities in (j4.5|) - (j4.7|) are 
Lipschitz continuous. We integrate the difference of the two equations ()6.ip . written for 
(0i,C/i,xi) and (^2,^^25X2)5 from to t, and test by w = Oi — 62- This yields 

f \9{x,t)\'^ dx + ^ (rq [ \VQ{x,t)\'^dx + Ri [ h{x){Q - erf{x,t) da{x) 

<^2(||^o|li2m)+ fmO?dC+ I h{x)\Q-Qr\\ev\{x,t)da{x) 
^ Jo Jan 




I an 

2 



n ^ Jo 



{mx,o\ + \ut{x,c)\ + mx,o\ + \uix,c)\)d^') dx). (6.2) 



Repeating the argument of [1U\ Proposition 4.3] or [H Proposition 3.4] about the L^- 
Lipschitz continuity of solution operators to gradient flows, we obtain for the solutions to 
(|i76|) - (|T71) for a.e. {x, t) E Qt the estimate 

[\\Xt{x,T)\ + \Ut{x,Tmx,Od^ + \x{x,t)\ + \U{x,t)\ (6.3) 
Jo 

<R3(^l i^o(e)|de + ixo(x)i + |{7o(x)i + ^ (i^(^,e)i + ^i^(y,oidy) d^). 

Integrating ()6.2p from to t and using (|6.3p together with Gronwall's argument, we 
obtain for each t € [0, T] the estimate 

|0(x,e)pdxde+ [ \x{x,t)\'^dx+ [ \U{x,t)\'^dx 
Jn Jn Jn 

< ^4(||^o|li2(f,) + WMUn) + I|f>o|li2(o) + ^ l^'o(6l' d^ 

+ t [ h{x)9l{x,0da{x)dC). (6.4) 
Jo Jan ' 

This concludes the proof of uniqueness of solutions and of Theorem W. 
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